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Abstract 

The Mathisson-Papapetrou-Dixon equations for a massive spinning test particle in plane gravitational 
waves are analysed and explicit solutions constructed in terms of solutions of certain linear ordinary differ- 
ential equations. For harmonic waves this system reduces to a single equation of Mathieu-Hill type. In this 
case spinning particles may exhibit parametric excitation by gravitational fields. For a spinning test particle 
scattered by a gravitational wave pulse, the final energy-momentum of the particle may be related to the width, 
height, polarisation of the wave and spin orientation of the particle. 

1 Introduction 

One of the most striking predictions of general relativity is the existence of gravitational waves. Such waves 
are thought to be produced by astrophysical phenomena ranging from the coalescence of orbiting binaries 
to violent events in the early Universe. Their detection would herald a new window for the observation of 
natural phenomena[l]. An essential mechanism by which gravitational waves interact with matter relies on 
the tidal forces produced on objects with structure in a gravitational field. This mechanism is responsible 
for forces and torques experienced by extended bodies or elementary particles with angular momentum in 
curved spacetime. By neglecting self-gravitation and back-reaction the dynamics of classical test particles 
with angular momentum was first studied in detail by Mathisson, Fock, Papapetrou et al [2]. The theory was 
further clarified by Dixon [3] using a rationalised multipole expansion technique and developed by Ehlers, 
Rudolph [4] and others. 

Up to the dipole approximation the motion of a classical spinning test particles is governed by the Mathisson- 
Papapetrou-Dixon (MPD) equations which encode a spin-curvature coupling that may play an important role 
in astrophysical contexts. In particular the MPD equations predict gravitational spin-spin interactions between 
rotating stars and orbiting massive spinning particles. The ever-increasing interest in black holes and the antic- 
ipation of gravitational wave astronomy have stimulated renewed interest in the dynamics of spinning particles 
near massive compact objects. This includes gravitational waves generated by spinning particles spiralling 
into black holes [6, 7, 8], chaotic behaviour of spinning particles in the Kerr metric [9, 10, 11] and numerical 
simulations for general orbital motions of spinning particles in stationary spacetime [12]. 

By contrast the dynamical response of spinning test particles to strong tidal forces produced by nearby 
gravitational events is less well understood although it had been suggested that the scattering of particles by 
plane gravitational waves provides a good local representation for such a physical process [13]. This problem 
has been recently considered from a number of different perspectives [14, 15, 16]. In this paper we offer a 
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new approach. We exploit the Killing symmetries of the plane gravitational wave metric and suitable space- 
time charts to construct a class of solutions in terms of certain linear ordinary differential equations. For 
harmonic gravitational waves this system reduces to a single equation of Mathieu-Hill type. By analogy with 
the parametric excitation of dynamical systems it is suggested that nonlinear spin-curvature interactions may 
exhibit parametric excitation by gravitational fields. In the case of a spinning test particle scattered by a square 
gravitational wave pulse, the final energy-momentum of the particle is shown to depend on the width, height, 
polarisation of the wave and spin orientation of the particle relative to a geodesic observer. 



2 Equations of Motion for Spinning Test Particles 

In the monopole-dipole approximation the MPD equations determine the world-line of a spinning test particle 
with tangent vector V, momentum vector P and spin 2-form s in terms of the spacetime metric g, its Levi- 
Civita coimection V and Riemannian curvature tensor R. ^ By introducing the metric duals p = P = g{P, —) 
and V = V, denoting interior operator (contraction with any vector V) on forms by ly [17] the MPD equations 
[3] can be written in a compact form as 

yyp = lyf (1) 



with the spin condition [5] 



and forcing term 



Vvs = 2p^v (2) 



zps = (3) 



f =^-^~\Rabh^s)e'' he"" (4) 



in terms of the exterior product A and Hodge map -k associated with a metric with the signature (—,+,+,+). 
Thus the "inverse Hodge map" denoted by [17] acts on any g-form uj according to 'k~^uj = — (— l)'?(^~'')*a; 
satisfying ★(★"^w) = -k~^{-kuo) = uj. The curvature 2-forms i?";, in any basis {Xg} with co-basis {e"} are 
related to the curvature tensor components R"'bcd by 2 R"'b{Xc, X^) = R"'bcd' (O) b,c,d = 0, 1, 2, 3). In terms 
of the components /„ = = e"(y), p" = e"(P), Sab = s{Xa, Xb) (4) becomes 



fa = - ^ RabcdV^ s'^'^ ■ 

It follows that [6] the velocity of the spinning test particle takes the following explicit form 

' P - T-7- X-w7 ■ (5) 



pf Pf \ 4.pfpc- Rhode s''" S'^ 

The freedom to normalise V ensures that p^ Vc is arbitrary and once a parameterisation of the world-line has 
been chosen (5) permits the computation of the world-line [4]. The norm of the momentum vector given by 



m = y/-giP,P) (6) 

may be identified as the mass of the particle. This is used to define the normaUsed momentum vector given by 

P 

U = - (7) 
m 



'In this paper we choose units where c = G = 1. 
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In Minkowski spacetime U coincides with the 4-velocity V/ ^/—giV-, V) but in general these two vectors differ 
in curved spacetime. Using the 1-form « = [/ the spin l-form I is defined by 

l = -\*{u^s) (8) 

and spin vector by L = 1. The angular-momentum 2-form s may be related back to the spin 1-form by 

s = 2iu-k''H. (9) 

The norm of the spin vector L denoted 



I = ^g{L,L) (10) 

defines the "spin" of the particle. For convenience we introduce the "reduced" spin vector , S = L/m, and 
reduced spin 1-form a = l/m such that the norm a = y/g{T,, S) = i/m denotes the spin to mass ratio. It 
follows from (1), (2) and (3) that a and m are constants of motion and the vectors U and S are orthogonal, i.e. 

5(C/,S) = 0. 

We are interested in solving the MPD equations in a gravitational wave spacetime that admits enough 
symmetry to permit the construction of constants along the world-line C. If sufficient constants can be found 
they enable one to integrate the equations of motion in terms of known functions. Suppose the metric g admits 
a Killing vector field K and let k = K.lt follows from the defining equation Ck g = that 

{VYk){Z) + {Vzk){Y) = (11) 

and 

Vy dfc = 2 A;" Rab (12) 
for any vectors Y and Z where fc" = e^{K) and y" = e'^(F). The preceding relation is equivalent to 

in component notation. By introducing the 4-form 

the MPD equations (1), (2) and (3) together with the above relations then imply that 



(13) 



1 1 

VvlK = Vvp/\*k+pA-kWv + - Vy (dfc) A*s + - d/c A *Vy s = 0. (14) 



Therefore for spacetime with a Killing vector K the quantity 

„ab 



(15) 



is preserved along the world-Une of a spinning test particle [18]. 
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3 The MPD Equations in a Plane Gravitational Wave Metric 

In the coordinate chart {t,x,y, z) we represent a plane gravitational wave metric g in the Kerr-Schild form 
[19]: 

g = H{jl, X, y) d/i (8> d/x — - (d/t (8) d^' + dz> (8) d/t) + dx dx + dy ® dy (16) 

where fi = t — z and v = i + z and these coordinates have dimension of length. We choose 

H{fi,x,y) = f{fi){x^-f) (17) 
for an arbitrary wave profile /(/t) and adopt the orthonormal co-basis 

1- H{il,x,y))^ 





_ di> 




~ T 




= dx 




= dy 




_ du 




~ T 



{1 + H{fi,x,y))^ (18) 

with dual basis {Xa} such that e"'{Xh) = S'^b. This basis is parallel along the geodesic observer O : t 
{i{r) = T,x = 0,y = 0, z = 0), i.e. S/ Xa\o = 0, and takes the form {Xq = d^, Xi = d^, X2 = dy, X^ = 
dz} which sets up a "local Lorentz frame" along O. In the coordinates {fi, u, x, y) this observer is parametrised 
by O : r 1-^ ipir) = r, z)(r) = r, f = 0, y = 0). 

It is instructive to study the motion of a spinning test particle with mass m, spin I = am and ini- 
tial location coinciding with the observer O at r = excited from rest by the incident gravitational wave. 
Thus we parametrise the world-line of this particle by C : A 1— > (/i(A) = A, i>(A),f (A),y(A)) for functions 
z>(A), f (A), y(A) subject to the initial conditions: 

i>(0) = x(0) = y(0) = (19) 
Let ii" = e"(f/) and a" = e"(E) be functions of A satisfying 

%6u"(A)u^(A) = -1 (20) 
VabCJ%\)a\X) = (21) 
%5^«(A)(7''(A) = (22) 
for 'qab = diag(— 1, 1,1,1). The initial conditions for are simply 

ii°(0) = 1, -ti^(O) = ^2(0) = ^3^0) = (23) 
while the initial conditions for cr" incorporate the initial spin orientation according to 

CT°(0) = 0, cr^(O) = ai, 0-2(0) = a2, a^(0) = (24) 
where constants a\,a2, 0:3 satisfy a\ + a\ + a\ = a^. 
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The Killing symmetry of g is more conveniently expressed in the coordinates {ji, v, x, y) defined by [19] 

ji = H (25) 

V = v + x'^ a{ii) a'{n) + 6(/i) (26) 

X = a{ii) X (27) 

y = b{fi)y (28) 
where the metric functions a{^) and 6(/x) satisfy 

a"(/i) = /(M)a(M) (29) 

b"{p.) = -f{ii)h{ii) (30) 
with a convenient choice of the (gauge fixing) conditions: 

a(0) = 6(0) = 1 (31) 

a'(0) = 6'(0) = 0. (32) 
In this coordinate chart the metric g takes the Rosen form 

g = (d^ (g) di/ + di/ ® d/i) + a(/i)^ dx (g) dx + 6(//)^ dy dy (33) 
admitting five Killing vectors [19]: 

d d d 

Ki = —, K2 = ^, Ks = — (34) 



du^ dx^ dy 

d ^ d 
ox ov 

d ^ d 

oy av 



^4 = ^(M)^ + 2a;|- (35) 



K5 = B{^) — + 2y— (36) 



where the functions A(iJ,), B{ij) satisfy 



For simpUcity we set A(0) = B(0) = 0. It follows from (15), (19), (23) and (24) that 

Tfb 

Ci = --, C2 = 0, C3 = 0, C4 = ma2, C5 = -mai. (38) 

From (19) and (25) - (28) in the coordinates {ji, u, x, y) the initial conditions for the world-line of the spinning 
test particle C : A {t^-V^) = ■> x{X) , y{\)) read 

m = x{0) = y{0) = (39) 

Thus (20) and (38) yield 

AA) = l + ^aW + f W (40) 

u^A) = -a2 a'(A) (41) 

u2(A) = ai 6'(A) (42) 

uHX) = fa'{Xf + ^b'{Xf (43) 



while (22) and (38) give 

- 2{y{\) - Qi)Q2i>(A)a'(A) - 2(i(A) + c.2)c.,o(A)6'(A)} (44) 
.>(A) = -(.(A) - cOKA) + l-2aW(xy(X) 

+2 02 a'(A) [(2/(A) - ai)a2KA)a'(A) + (x(A) + a2)aia(A)6'(A)] } (45) 

.^(A) = (.(A) + «2)a(A) + ,.,.(,)2^ ^'(A)c^^(A) 

-2 ai 6'(A) [(y(A) - ai)a2biX)a' (X) + (x(A) + a2)aia(A)5'(A)] } . (46) 
Furthermore (21) and (38) imply quadratic equation for (7^(A): 

a^{Xf + V{X) u^{X) + Q(A) = (47) 

where 

V{X) = -2 {y{X) - ai)a2b{X)a\X) - 2 (x(A) + a2)aia(A)6'(A) (48) 



and 



Q{X) = -^{ala'{Xf + aly{Xf-2Y 

^ {y{X) - aifb{Xf s^2alala'{Xfb'{Xf + a^a'iXf + a\b'{Xf - Aa\b'{Xf + 4} 



(a;(A) + a2)MAf |2a2a2^'(A)26'(A)2 + aia'(A)4 + a\b'{Xf - Aala'{Xf + 4} 



+2aia2o(A)6(A)a'(A)6'(A)(x(A) + a2)(y(A) - ai). (49) 

The differential equations for z/(A), x(A), y{X) are obtained from substituting (7), (9) and (40) - (46) into (5) 
as follows: 

^x(A) = 2 /(A) {aia(A)6(A)6'(A)(x(A) + «2)(y(A) - ai) + a26(A)'a'(A)(y(A) - aif 
-{y{X) - a,)b{X)aHX)} / {a{X) [al a'{Xf + a\ b'{Xf - 2) } - (^(A) + «2)a^(A) ^^^^ 

Ay(A) = 2 /(A) {a2a(A)6(A)a'(A)(x(A) + a2){y{X) - «i) + aia(A)26'(A)(a;(A) + as)' 
-(x(A) + a2)a(A)a-^(A)} / {6(A) {al a' {Xf + a? ^'(A)^ - 2) } - ^^^^^ 7 "'^^'^^^ (51) 



Ai.(A) = 1 + (x(A) + a-2fa'{\f + (y(A) - aifh'{\f 



+/(A) {Aa\\){x{\) + a2){y{\) - ai)(a(A)6'(A) + 6(A)a'(A)) 



+2 (x(A) + a2fa{\f (2 - al a'{Xf + a? 6'(A)2 - 2ai 6'(A)\(A)) 
-2 (y(A) - ai)%{\f (2 + ai a' {\f - a\ h'{Xf + 2^2 a'(A)\(A)) 
-4a(A)6(A)a'(A)6'(A)(x(A) + a2)(y(A) - ai){aix{\) + a2y(A))} 

/(a^a'(A)2 + a2 h'{\f-2). 



(52) 



Note that i/(A) is obtained by directly integrating the R.H.S of (52) which is in turn decoupled from (50) and 
(51). Thus the world-line of the spinning particle is effectively determined by solving two coupled non-linear 
differential equations (50) and (51). 

Given the initial reduced spin vector components a\, 02, 03 with respect to the orthonormal basis {Xg} 
along C the evolution of the normalised momentum vector U is readily available from the explicit expressions 
(40) - (43) for some gravitational wave profile /(/x) and the corresponding a(/x) and (cf. (29), (30), (31), 
(32)). However to obtain the world-line C the differential equations (50), (51), (52) must be solved together 
with (47) subject to initial conditions (39) and the evolution of the reduced spin vector S follows from (44) - 



A class of trivial solutions is obtained if the spin vector is initially in the direction of the incident gravitational 
wave (tti = Q!2 = 0, q; = |q;3|). In this case the world-line C stays coincident with the geodesic observer 
O while both U and E remain covariantly constant along the world-line. A class of non-trivial solutions 
describing non-geodesic motion of a particle with a spin vector which stays parallel along its world-line and is 
transverse to the propagation direction of the gravitational wave may be obtained as follows: Suppose a = |ai |, 
Q!2 = "3 = and 



(47). 



4 Solutions with Parallel Spin Vectors 



x(A) = 0, (j3(A) = 



(53) 



then (50) is immediately satisfied and (51) becomes 



6(A) 



(54) 



which has a solution 




(55) 



compatible with the initial condition 2/(0) = 0. Equation (52) then reduces to 



where (30) has been used. With the initial condition v{Q) = this can be integrated to yield 




(56) 




(57) 
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where 

= tfrnv- (58) 
Jo 

Substituting the above equations into (40) - (46) we have 

u^{\) = 1 + ^b'{\f, u\\) = 0, u\\) = ai b'{\), u^{X) = ^ b'iXf (59) 

cj°(A) = 0, cj^(A) = ai, a^{X) = 0, cj^(A) = 0. (60) 
Using (25) - (28) solutions (53), (55), (57) in the coordinates {i, x, y, z) are given by 

i(A) = A + ^{(6(A)-2)6'(A) + F(A)} (61) 

x(A) = (62) 
y{\) = ai (6(A) - 1) (63) 

^(A) = f {(6(A)-2)6'(A)+F(A)}. (64) 

A further class of solutions describing non-geodesic motion is associated with the choice ai = as = 
a = |q;2 I and may be derived in a similar way: 

2 2 

uO(A) = 1 + ^a'(A)^ u\X) = -a2 a'(A), ^^(A) = 0, u\X) = ^ a'{Xf (65) 

(T°(A) = 0, (T^(A) = 0, (7^ (A) = a2, (T^(A) = (66) 
with the corresponding world-line given by 

iW = A + f {(a(A) - 2) a' (A) - F(A)} (67) 
x{X) = -a2 (a(A) - 1) (68) 
y{X) = (69) 



^(A) = ^{(a(A)-2)a'(A)-F(A)} (70) 

in the coordinates {i,x,y, z). 

The above solutions are valid for arbitrary wave profile /(/x) and the corresponding functions o(//), 6(/u), 
F(/x) satisfying (29), (30), (31), (32) and (58). For example a harmonic gravitational wave with angular 
frequency u, dimensionless amplitude h and a mean value k/2 may be represented by 

/(//) =a;2 1^-^ cos(6c;m)1. (71) 



^4 2 

Introducing the variable ( = uj iJ./2 one immediately sees that (29), (30) can be cast into the canonical form of 
the Mathieu equations [20]: 

^"(t;) ={^-2^cos(2C)}a(^^) (72) 



6 { — ) = {-fe + 2/icos(2C)} 5 ( — ) . (73) 



Thus a{fi) and b{fi) may be expressed in terms of Mathieu functions with characteristic numbers ±k and mod- 
uli ±h respectively. This suggests that nonUnear spin-curvature interactions may exhibit parametric excitation 
by gravitational fields [21] in analogy with the parametric excitation of dynamical systems. 



5 Scattering by a Gravitational Wave Pulse 



We now consider a spinning test particle scattered by a square gravitational wave pulse represented by the 
metric (33) with metric functions a(/x) and associated with the wave profile /(//): 



if ^ < 
l/C"^ ifO<n<S 
otherwise. 



(74) 



Since both the width 5 and height parameter C carry dimension of length it is convenient to choose a length 
unit in which C is unity. To recover from this non-dimensionalisation one simply applies the substitutions: 
X X/jO, II ^ li/C,a ^ a/jC,m m/£, etc. With £ = lit follows from (29), (30), (31), (32) and (74) 
that 



1 if ^ < 

cosh(^) if < /i < (5 

, sinh((5) jjL — sinh((^) d + cosh((5) otherwise 



(75) 



1 if ^ < 

cos(^) if < /i < (76) 

■ sin(5) /J, + sm{6) 6 + cos{6) otherwise 



As is well known such a sandwich wave gives rise to two half Minkowski spacetimes namely for ^ < 
and A4-^- ioi fj, > 5 corresponding to the domains before and after the passage of the gravitational wave pulse. 
Within each half Minkowski spacetime the orthonormal basis {Xq = df,Xi = dx,X2 = dy,X^ = dz} 
defines a "global" Lorentz frame, whereas {^o} along the geodesic observer O still remains a local Lorentz 
basis in the "sandwiched" domain V where < /x = t — i < (5. It is therefore interesting to analyse how a 
spinning test particle initially following O with a spin orientation with respect to {Xg} is scattered by the wave 
pulse and establish the final U and S. In M.- and M.+ these vectors have different constant components u"', 
cr" and changes in these quantities along C across the domain V represent the gain in energy-momentum and 
precession of spin in the Lorentz frame associated with the observer O. To establish the relations describing 
these changes consider the evolution of the components of U which is readily evaluated from (40) - (43) to be 



Domain: 


M- (A < 0) 


V{0<X<S) 


M+ (A > S) 


vP{X) 


1 


l + ^sinh2(A) + ^sin2(A) 




u\X) 





— a2 sinh(A) 


—a-i sinh((5) 


u\X) 





— cti sin(A) 


—a\ sin((5) 


u^X) 





f sinh^(A) + f sin^(A) 


Msinh2(<5) + ^sin2(5) 



(77) 
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using (75) and (76) without explicitly solving (50) - (52). Since in Minkowski spacetime the normalised vector 
U coincides with the particle's 4- velocity, for A > 5 the motion of the particle in A^_|_ has a 4- velocity with 
constant components m'*(A) = u"'{S) given in (77). Thus the difference between these values and the initial 
components = l,u^ = = = determines the deflection of the particle motion relative to the observer 
O. 

To obtain the evolution of the reduced spin vector it is necessary to solve (50) - (52). If the initial reduced 
spin vector is perpendicular to the propagation direction of the wave then as described in the previous section 
these equations yield solutions (59) - (64) and (65) - (70) where a(/i), b{fi) are now fixed by (75), (76). For 
general ai,a2, as, a power series solution of (50) - (52) for < A < 5 may be obtained and from (44) - (47), 
(25) - (28) together with matching the particle trajectory C across domains M-yV, M.+ this yields 



Domain: 


M_ (A < 0) V{0<X<S) 


M+ (A > S) 




£(A) 


A 


(l + a?-a^)A + 0(A=^) 


t{S)+u^{d) (A- 


S) 


x{\) 





-ai aa A + Y (a^ - 3a? - ai - 1) A^ + 0{X^) 


x{d)+u^{5) (A- 


6) 


m 





a2 as A - y (a^ - 3ai - a? - 1) A^ + 0{X^) 


y(,5) + «2(<5) (A- 


6) 


i(A) 





{al - ai) A + O(A^) 


z{S) + u^i6) (A - 





(78) 



and 

Domain: M- (A < 0) T> {0 < X < d) M+ (A > 6) 



aO(A) 





-2 ai a2 A + ^ (a? + ai) A^ + O(A^) 


a'{6) 


aHX) 


ai 


a, + ^{a^-al + a^^)X' + OiX^) 


a\S) 


a\X) 


a2 


a2 + ^{a' + al-al)X' + 0{X') 




aHX) 


as 


as - ^ (a? + a|) A^ + O(X^) 





(79) 



As expected for A > (5 the "out-going" motion of the particle in M+ follows a geodesic with constant 4- 
velocity components u^{6) given in (77) and constant reduced spin vector components (t"'{S) (given approxi- 
mately in (79)). This represents the precession of the spinning particle due to the scattering by the gravitational 
wave pulse. 

6 Conclusions 

We have considered the effects of a plane gravitational wave spacetime on the motion of a spinning test particle. 
A class of solutions of the MPD equations has been constructed in terms of solutions of certain linear ordinary 
differential equations. The motion of the particle with a parallel spin vector transverse to the direction of 
propagation of a polarised harmonic gravitational wave has been formulated in terms of Mathieu functions. 
By analogy with the parametric excitation of dynamical systems it is suggested that nonlinear spin-curvature 
interactions may exhibit parametric excitation by gravitational fields. For a spinning test particle scattered by 
a gravitational wave pulse the out-going energy-momentum of the particle has been expressed in terms of the 
width, height, polarisation of the wave and spin orientation of the particle relative to a geodesic observer. 
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